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In spite of the apparent similarity of micro-branching instabilities in different brittle materials, 
we propose that the physics determining the typical length- and time-scales characterizing the post- 
instability patterns differ greatly from material to material. We offer a scaling theory connecting the 
pattern characteristics to material properties (like molecular weight) in brittle plastics like PMMA, 
and stress the fundamental differences with patterns in glass which are crucially influenced by 3- 
dimensional dynamics. In both cases the present ab-initio theoretical models are still too far from 
reality, disregarding some fundamental physics of the phenomena. 



More than a decade of relatively high precision experi- 
ments on rapid crack propagation in brittle materials has 
ignited the interest of the physics community [l[ . It ap- 
pears that one important prediction of "Linear Elasticity 
Fracture Mechanics" 2], stating that cracks will accel- 
erate smoothly as they lengthen until they reach their 
asymptotic velocity (bounded either by the loading con- 
ditions or by the Rayleigh speed cr), is in fundamental 
contradiction with experiments. In thin brittle plates of 
width W (see the typical geometry in Fig. [I] W is in the 
z direction, with length in x and height in y much larger 
than W), at some finite fraction of the Rayleigh speed, 
the smooth dynamics of the crack evolving in the x di- 
rection is marred by the appearance of micro-branches in 
the x — y plane that do not cut through the whole width 
W. These micro-branches have a typical length £5, a 
typical width Az < W and they appear to recur with 
an average periodicity. We will refer to this recurrence 
scale as "the scale of noisy periodicity" Ax. The nature 
of this length-scale and the physics that determine it are 
the subjects of this Letter. In glass this scale appears to 
range between the measurement resolution (~1 fim) and 
a few mm 4], while in PMMA between tens and hun- 
dreds /im |3(, depending on the material properties and 
the experimental conditions. 

As is well known, linear elasticity fracture mechanics 
has no typical scale, and thus the appearance of such a 
scale is yet another demonstration of the need to modify 
this theory in the context of dynamic fracture. Indeed, 
there exist concepts in the classical theory of crack prop- 
agation that can be used to form a length scale. Recall 
that in the classical theory one asserts that as a crack of 
length L propagates at velocity v under loading condi- 
tions a, the energy G(a,L,v) released from the stressed 
material streams into the tip region where it is compen- 
sated by the total dissipation T(v), 



G(a,L,v)=T(v) 



(1) 



The dissipation function T(v) is not computable from 
elasticity theory, but it can be measured in experiments, 
(essentially by using Eq. QJ). Equipped with this func- 
tion and the Young modulus E of the material we can 




FIG. 1: Typical geometry of a dynamically generated crack, 
with the side branches that result from the instability leaving 
their mark on the faces of the crack. 



form a (velocity-dependent) scale from 
£(v) = T(v)/E . 



(2) 



To get a rough estimate of the resulting scale, we use 
the data for PMMA determined in |l| near the onset of 
the microbranching instability: T(v c ) ~ 3000 J/m 2 and 
E ~ 3 • 10 9 N/m 2 . This predicts £(v) in the fim range; a 
similar estimation for glass yields a nm scale. We thus 
understand that the scale computed in Eq. (0) may be 
relevant for PMMA, but appears utterly out of range 
for glass, where the observed scales of noisy periodicity 
must be sought elsewhere, as we discuss below. Given 
information about the largest stress value that marks the 
breakdown of the linear theory, say <r m , another length- 
scale can be constructed as 



£(v) = ET{v)/a 2 n 



(3) 



We estimate <j m by the yield stress of the material. 
For the materials under discussion the resulting length- 
scale is typically two order of magnitude larger than @ , 
putting it right in the range of scales of noisy periodic- 
ity for PMMA, but still off range for glass. We continue 
now to assess the relevance of (0 or (PJ to PMMA. Obvi- 
ously, the dissipation T(v) in PMMA should be a strong 



2 



function of the molecular weight M of the polymers that 
make the material, due to the increased density of entan- 
glements of the polymer chains. Indeed, a calculation of 
the M dependence of T in quasi-static conditions exists 
H, with the resulting prediction 



T(v = 0)~e 



-Mq/M 



(4) 



where Mo is the smallest molecular weight below which 
the polymers do not contribute to T. On the other hand, 
the Young modulus E and the yield stress were shown 
0,0 to be M-independent. 

The length (or [3J) will be dressed by dynamical 
effects, which we introduce in the form 



i(v) = e- M °' M f{v/c R ) 



(5) 



To proceed, we adopt the scaling assumption that £(v) is 
the only typical scale in the problem, characterizing the 
length £b of micro-branches in the x — y plane as well as 
Ax in the x — z plane. In addition, we will assert that 
in PMMA (in contrast to glass, and see below for fur- 
ther discussion) the relevant length-scales are determined 
by a competition between the side micro-branch and the 
main crack, without important dynamical coupling to the 
third dimension. In this respect the mechanism underly- 
ing the repetitive nature of the microbranching process 
in PMMA is qualitatively similar to the 2-dimensional 
branching model developed recently in 8]. In this model, 
the local competition between the side branch and the 
main crack temporarily slows the main crack below the 
relevant branching velocity. When the main crack regains 
the critical velocity another side branching event occurs 
(see Fig. 7 in 8]). It was found (see Fig. 6 in 8]) that 
the velocity dependence of has the following form 



h(v) = £ b (v c ) [1 + ol{v - v c )/c R ] 



(6) 



Here £b(v c ) is the finite length of the side branch which 

occurs at threshold and a is a dimensionless coefficient. 

i—i 

In |2| this length was introduced by hand to represent 
the scale over which linear elasticity theory is not appli- 
cable (the so-called process zone); here we take it as the 
undressed length proportional to exp {—Mq/M). Collect- 
ing our assumptions together we present the prediction 
for the scale of noisy periodicity and the length of the 
micro-branches in PMMA: 



Ax(v : M) ~ £ b (v, M) ~ e - M °/ M 



1 



OL- 



CR 



• (7) 



Note that if indeed the scale £b(v c ) is physically identified 
with the process zone, it may have a velocity dependence 
of its own. If this velocity dependence is linear, the pre- 
diction remains unchanged. 

It is possible to test this major prediction against avail- 
able experiments. The influence of the mean molecular 
weight M on the scale of noisy periodicity in PMMA was 



studied in at a given average velocity. Fig. [3 displays 
their experimental data for M log (Ax) as a function of 
M (taken from Table 1 in 9]). The almost perfect lin- 
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FIG. 2: M log (Ax) as a function of M. The filled circles are 
the experimental data taken from |9| (see Table 1) and the 
solid line is a linear fit. 

earity of the graph lends a strong support to the M de- 
pendence in our scaling law. Next, the branch length at 
a fixed mean molecular weight as a function of the nor- 
malized mean velocity for PMMA is shown in Fig. 03 It 
is clear that £ b is linear in v as predicted by our scaling 
relation. Moreover, the fit indicates the existence of a fi- 
nite length £(v c ) « 30/im as predicted by the scaling law. 
This finite length is determined uniquely in this graph 
once v c is known. The value of £(v c ) is confirmed in 0. 

In summary, the scale of noisy periodicity in PMMA 
appears to be adequately described by a scaling theory 
in which the main scale depends on the material proper- 
ties T, E and possibly <r m , dressed by dynamical effects 
which are not fundamentally dependent on coupling to 




FIG. 3: The branch length £b as a function of the mean nor- 
malized velocity for PMMA. Data taken from 10]. The solid 
line represent the liner fit with a ~ 60, v c ~ 0.365cr, and 
£(v c ) « 30/xm. 
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dynamics in the third dimension. There is nothing like 
that in glass. First, we remarked before that the length- 
scales © and © are in 1-100 nm range, whereas the 
noisy periodicity in glass appears on scales up to a few 
mm. Moreover, we will point out now that there is no 
typical scale of noisy periodicity in glass experiments; 
the distances between micro branching events can vary 
over three orders of magnitude (from the measurement 
resolution to the scale of the width W) and they are de- 
termined by a dynamical mechanism that rests crucially 
on the dynamics of the crack front in the narrow third 
dimension along the width W . 

In glass, once a micro-branch forms at a location z 
along the width of the sample and with a width Az, the 
next branching events fall very precisely along a strip of 
width Az localized at the same location z; see Fig. 0]as 
an example. The crucial phenomenon that is responsible 
for these highly correlated structures is the interaction 
of the crack front with spatial heterogeneities (regions 
of variable T). These can generate non-decaying waves 
that travel with velocity cpw (relative to the heterogene- 
ity) on the front as it propagates 11]. Plane numerical 
simulations have shown that the interaction of the 
front with an asperity leads to a velocity overshoot ex- 
actly ahead of the interaction site, see Fig. 



the order of the width of the micro-branch Az. Since the 
velocity of the front wave along the front is \/ c 2 FW — v 2 
one obtains 



Az(v) 



(8) 



FW 



Therefore, the noisy periodicity scale Ax is estimated as 

Az(v) v 



Ax(v) 



(9) 



FW 



Since the combination vj \/ c 2 FW — v 2 is a slowly in- 
creasing function of v (0.94cr<cfw(v)<cr) in the 
relevant range of velocities, one expects the ratio 
Ax(v)/Az(v) to be nearly constant and this is indeed 
the case as shown in |4|, where the relation is stated to 
be Ax ^4-5 Az. The point that we want to stress here is 
that in contradistinction with PMMA, in glass there exist 
no typical scale for the noisy periodicity. Ax(v) can be 
anything, depending on the local asperity size that nu- 
cleated the first micro branch of size Az(v). After that 
the dynamics will sustain repeated micro-branches every 
Ax c^A-bAz. The only thing that can happen is that 
Az(v) might increase as a function of v, which is the last 
issue that we discuss. 





FIG. 4: The correlated micro-branches in glass samples. Each 
series of micro-branches is localized in a strip of width Az, 
and has a noisy periodicity Ax~4-5Az. Adapted from 



Since the creation of a micro-branch increases T lo- 
cally, one expects the generation of front waves, as in- 
deed was observed in glass and can be seen in the inset 
of Fig. EI The resulting velocity overshoot just ahead 
of the micro-branching event was interpreted by Sharon, 
Cohen and Fineberg U as the source of the well-defined 
lines of micro-branches in glass. In this picture, the cre- 
ation of a micro-branch reduces locally the velocity of 
the main crack below the micro-branching velocity, but 
this velocity is reached again due to the velocity over- 
shoot just ahead. The typical time r for this process is 
the time needed for the front wave to travel a distance of 




FIG. 5: The front waves generated by the interaction of the 
front with an asperity. The front velocity profile is plotted 
at constant time intervals. Note the velocity overshooting 
just ahead of the asperity. d,s,R are the cones generated by 
the dilatational, shear and Rayleigh waves, respectively. The 
relation cfw ~ cr is apparent. Adapted from [lg. 

As noted earlier, Eq. [I] predicts that cracks acceler- 
ate smoothly towards their asymptotic velocity. A cru- 
cial assumption leading to this prediction is that a single 
crack whose front is essentially a point (there is trans- 
lational invariance in the z-axis) can expend the extra 
energy flowing to its tip by adjusting its velocity and 
increasing the kinetic energy of its surroundings. This 
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assumption holds up to v c , but breaks down for v > v c . 
At the onset of instability a multiple crack state is formed 
by repetitive frustrated micro-branching events and these 
micro-branches break the translational invariance in the 
third dimension, occupying a finite width Az < W. The 
essence of the instability is that the additional energy 
supply is expended on creating more surface area per unit 
crack extension, even at a constant velocity. This pro- 
cess introduces new dynamic variables to the problem: 
the number of coexisting branches n, the morphology of 
the branches in the x-y planes and their width Az. All 
these variables contribute to the additional energy dissi- 
pation per unit crack extension. It was shown that 
the total surface energy created per unit crack extension 
is proportional to the total energy release rate, implying 
that in this regime the fracture energy T(v) is constant, 
independent of the velocity. Assuming that Az charac- 
terizes the width of the micro-branches as long as they 
exist, we can write 
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FIG. 6: The pattern width Az as a function of the normal- 
ized mean velocity v/cr for PMMA (filled symbols) and glass 
(empty symbols). Adapted from |T{J. 



WG(v) = Az(v)f(n(v), ...)T + WT , 



(10) 



where the right hand side is a sum of the dissipative con- 
tributions from the micro-branches as well as from the 
main crack, respectively. The function f(n(v), ...) de- 
pends on the number of coexisting micro-branches, their 
morphology and maybe other dynamic variables. The 
last relation can be rewritten as 



Az(v) = f(n(v),...)W 



G(v) 



1 



(11) 



Although f(n(v), ...) is not known in adequate detail, 
making the last equation somewhat tentative, we can 
still draw relevant qualitative conclusion. Since G(v) is 
known from experiments to be a strongly increasing func- 
tion of v 13], we can safely expect that Az(v) is also a 
strongly increasing function of v. This expectation is sup- 
ported by the data; the dependence of Az on the mean 
velocity v for both PMMA and glass is shown in Fig. E3 
We stress that in glass the increase in Az(v) must affect 
the scale of noisy periodicity, whereas in PMMA it does 
not. In glass, the dynamics along the third dimension are 
strongly coupled to the dynamics in the longitudinal di- 
rection (by the front waves mechanism) and the dynamic 
length-scale Az controls all the other length-scales. In 
PMMA, the dynamics along the third dimension are only 
weakly coupled to the dynamics in the longitudinal direc- 
tion, where the length-scales are determined by material 
properties and the elastodynamic competition between 
the micro-branches and the main crack. 

In summary, models of dynamic fracture that are de- 
veloped based on the continuum equations supplemented 
by criteria for crack evolution can hope to adequately de- 
scribe the dynamic instabilities only if the physics that 
determines the fundamental scales were incorporated. In 
glass it appears that without adequate treatment of the 



third dimension and the physics of the crack front, not 
much insight into the micro-branching process can be 
gained. In PMMA the coupling to the third dimension 
might be weaker, but one should resolve how the ma- 
terial parameters appear in the continuum description. 
Completely different scaling laws appear in these two ma- 
terials due to the different mechanisms that control the 
scale of noisy periodicity; probably other materials may 
increase the richness of the phenomenology. It might be 
worthwhile to stress less the presence of universal features 
and more the interesting physics that can be learned even 
with simple scaling theories. 
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